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The structure and vibrational spectrum of Gd3N@Cso is studied through Raman and inelastic 
electron tunneling spectroscopy (IETS) as well as density functional theory (DFT) and universal 
force field (UFF) calculations. Hindered rotations, shown by both theory and experiment, indicate 
the formation of a GdsN— Cso bond which reduces the ideal icosahedral symmetry of the Cso cage. 
The vibrational modes involving the movement of the encapsulated species are a fingerprint of 
the interaction between the fullerene cage and the core complex. We present Raman data for the 
Gd 3 N@C 2 n (40 < n < 44) family as well as Y 3 N@C 8 o, Lu 3 N@C 8 o, and Y 3 N@C 8 8 for comparison. 
Conductance measurements have been performed on Gd 3 N@Cso and reveal a Kondo effect similar 
to that observed in Ceo- 
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I. INTRODUCTION 

Endohedral fullerenes (endofullerenes), due to their 
unique electronic and geometric structure, have been at- 
tracting the interest of chemists, physicists, and mate- 
rial scientists since they were first synthesized in 1985 
Fullerenes can act as hosts that encapsulate other 
atoms, molecules, or atomic clusters. The resulting 
endofullerenes have initiated considerable research be- 
cause the entrapped atoms bring with them an array 
of useful physical properties. The potential applica- 
tions include electronic devices Q, organic solar cells 
3[, spin-based quantum computing \4, l5| , and medicine 
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specifically contrast agents for magnetic resonance imag- 
ing (MRI) as well as radiopharmaceutical candidates 
@, @. Endofullerenes created by an arc- vaporization 
technique [22j] have been known to encapsulate up to 
three individual metal atoms as well as clusters of the 
form M 3 N, M 3 C 2 , M 2 C 2 , and M 4 2 (M = Sc, Y, and 
several lanthanides) The possible isolation of mag- 
netic or radioactive molecules from the surrounding envi- 
ronment have made endofullerenes interesting for many 
fields. In particular, the Gd 3 N@C 2n family has received 
much interest due to their application as MRI contrast 
agents and the confinement of three paramagnetic Gd 
ions in a relatively small space (~ 1 nm). 

The vibrational modes involving the movement of the 
encapsulated species are a fingerprint of the interaction 
between the fullerene cage and the metal complex. Based 
on the mode frequencies, the experimental value of the 
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metal— cage bond strength can be determined and the 
amount of char ge t ransfer from the cage to the core can 
be investigated [23l H3, [H| . It has been shown that some 
isolated fullerene cages, though electronically unstable, 
can bind strongly to the metal core and become stable 
(26[. Although certain gadolinium endofullerenes have 
been isolated and studied by X-ray diffraction, vibra- 
tional spectroscopy has not been performed on a number 
of species Q. 

In this paper, we present a detailed Raman and 
inelastic electron tunneling spectroscopy (IETS) anal- 
ysis of Gd 3 N@C8o, as well as conductance measure- 
ments demonstrating a Kondo effect. Raman analysis 
of Gd 3 N@C 2n (40 < n < 44) has been performed, as well 
as Y 3 N@C 8 o, Lu 3 N@C 8 o, and Y 3 N@C 8 s for compari- 
son. Vibrational assignments of the Gd 3 N@C8o molecule, 
Ih — C®Q cage, and intrinsic Gd 3 N cluster were carried out 
by density functional theory (DFT) and universal force 
field (UFF) calculations. Hindered rotations, shown by 
both theory and experiment, indicate the formation of 
a Gd 3 N— Cso bond which reduces the ideal icosahedral 
symmetry of the Cso cage. 



II. BACKGROUND 

An isolated Gd atom has a large magnetic moment 
of 8 where 7 hb is derived from localized / elec- 
trons [26[. When embedded inside a fullerene, the Gd 
atom bonds through the delocalized s and d electrons 
and the localized / states maintain their atomic spin 
moment [26j|. The valence state of gadolinium has been 
determined experimentally by ESR (Electron Spin Res- 
onance), XPS (X-ray Photoemission Spectroscopy) and 
X-ray diffraction to be Gd 3+ (27|, [H, S El [M 32] . For 
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the Gd3N@C2n family, Gd can formally only be in the 
+ 1 charge state, because with three additional electrons, 
one from each Gd, the nitrogen is saturated in a formal 
—3 charge state. However, the highest charge state is 
Gd 3+ and therefore the whole GdsN cluster can donate 
six electrons to the cage resulting in (Gd 2+ )sN(Q)C2~. 

In the case of endohedral doping, the stability and sym- 
metry of the cage depends on the properties of the encap- 
sulated species, mostly on the amount of charge transfer. 
As a consequence, the isomer favored upon incorpora- 
tion of an atom is usually different than the stable iso- 
mer of a hollow cage. Considering GdsN@C8o, the 1^ 
isomer of an isolated Cso cage is unstable, having a four- 
fold degenerate HOMO occupied by only two electrons. 
The empty cage would undergo a Jahn- Teller distortion, 
leading to an occupied non-degenerate HOMO and a low- 
lying triply degenerate LUMO. However, the transfer of 
six additional electrons from the three Gd 3+ atoms stabi- 
lizes the Ifc—Cso by completely filling the HOMO orbital 
[33( | . It is interesting to note that the intrinsic GdsN 
cluster and the Cso cage do not exist independently, but 
when combined form a stable molecule with a closed shell 
structure. 

The structure of GdsN® Cso has been determined by 
X-ray diffraction and the GdsN unit is found to be pyra- 
midal within the Ih cage, with the three Gd atoms posi- 
tioned over the centers of hexagons (34|. Unlike light 
mass pyramidal molecules, such as NH3, the N atom 
is not able to tunnel through the Gd3 plane due to a 
high energy barrier. However, when the intrinsic GdsN 
molecule is encapsulated in the Cso cage the charge cou- 
pling effectively lowers the energy barrier to ~ 91 meV, 
which allows tunneling to be possible [26[. As the clus- 
ter has a permanent electric dipole moment, the wagging 
mode involves a fluctuation in the overall dipole moment, 
which can be observed in the Raman spectra 85 cm -1 ) 
0. At room temperature, GdsN^Cgo is expected to be 
paramagnetic with the spin fluctuating between two dif- 
ferent multiplicities (M = and 22) [341 ]. 



III. EXPERIMENT AND RESULTS 
A. Spectrophotometry and Raman Measurements 

UV-Vis spectrophotometry data for Gd 3 N@C 2 n (40 < 
n < 44) are shown in Fig. 1. The samples were iso- 
lated by HPLC and the purity was confirmed by mass 
spectrometry. The endofullerenes were suspended in 
toluene (99.999%, Sigma-Aldrich[55|), which was back- 
ground subtracted, and plotted as a function of wave- 
length. 

The spectral onset for GdsN@C8o is around 940 nm. 
The HOMO-LUMO transition has a doublet structure 
with absorption maxima at 680 and 710 nm. The 
strongest visible absorption of GdsN^Cso is at 415 nm, 
and a shoulder is apparent at 555 nm [35]. The spectral 
onset for GdsN@C82 is around 950 nm and there are ab- 
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FIG. 1: Spectrophotometry data of Gd 3 N@C 2 n (40 < n < 44) 
in toluene. Unusually high absorption peak in Gd3N@Cs6 is 
observed. 



sorption peaks at 630 nm and 505 nm. The spectral onset 
for GdsN@C84 is around 985 nm and there are absorp- 
tion peaks at 915 nm and 635 nm. The spectral onset for 
GdsN@C86 is around 1005 nm and there is an unusually 
strong absorption peak at 640 nm, as well as one at 430 
nm. The spectral onset for GdsN@Cs8 is around 960 nm 
and there are absorption peaks at 740 nm, 545 nm, 460 
nm, and 405 nm. 

For spectroscopic studies, 10 fig of GdsN@C2 n (40 < n 
< 44) were suspended in carbon disulfide (CS2, 99.999%, 
Sigma- Aldrich) and used to dropcoat a gold-covered sil- 
icon substrate. The resulting polycrystalline films were 
dried under ambient conditions and placed in a Linkam 
liquid nitrogen cryostat stage. Raman spectra were stud- 
ied with 632.8 nm excitation from a HeNe laser. The 
scattered light was collected in the backscattering geom- 
etry by a triple-axis spectrometer (Jobin Yvon Horiba, 
T64000) equipped with a CCD detector. A resolution 
of 0.7 cm -1 (1800 gr/mm) was used for all Raman mea- 
surements and all measurements used 2.5 mW/cm 2 laser 
power, 30 min. accumulation time, and were performed 
at 90 K. 

The Raman spectra observed for GdsN^Cgo can be 
roughly divided into four parts: tangential Cso modes 
are found between 1000 and 1600 cm -1 , a gap-like region 
from 815 to 1000 cm -1 , radial Cso modes between 200 
and 815 cm -1 , and a fourth group from 200 cm -1 and 
below which has a counterpart only in the spectra of 
endohedral and polymeric fullerenes which exhibit low 
energetic Gd— cage, inter molecular, and center of mass 
(COM) modes. Low-energy Raman modes of GdsN^Cgo 
have been observed confirming the coupling between the 
core and cage by hindered rotation modes (Fig. 2). 

Several prominent modes can be identified such as the 
Cso H p (l) squashing mode at 234 cm -1 , the H^(2) twist 
mode at 361 cm -1 , and the A^(l) breathing mode at 430 
cm -1 . Below 200 cm -1 , several peaks are observed indi- 
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FIG. 2: Analysis of low-energy Raman lines of Gd3N@Cso 
taken at 90 K indicating a hindered rotation due to the 
coupling of the core complex to the cage. Inset gives the 
Gd 3 N@C 8 o structure. 



eating a hindered rotation due to the bonding of GdsN 
to the C8o cage. Other groups have seen evidence for 
the formation of a M3N— Cso bond which induces a sig- 
nificant reduction of the ideal 1^— Cgo symmetry [36|. A 
peak at 83.6 cm -1 appears to be the wagging mode of 
the nitrogen atom through the Gd3 plane, predicted by 
theory [26|. Other COM modes are identified, where the 
cage— core complex translates relative to the other. The 
frequencies of these peaks were compared with the other 
gadolinium endofullerenes and are listed in Table I. 

A similar procedure was performed to analyze the vi- 
brational spectrum of the Gd 3 N@C 2 n (41 < n < 44) 
family. It can be clearly seen that the low-energy Ra- 
man modes are changed as the number of cage atoms 
increases (Fig. 3). A common mode for all samples aris- 
ing from the coupling of the motion of the Gd complex 
and the fullerene cage seems to be apparent. The mode, 
defined from 140—165 cm -1 and identified boldly in Ta- 
ble I, can be assigned to the Gd— cage stretching mode. 
This particular mode is similar to the vibrational mode 
reported for Gd@C82, which exhibits a Raman frequency 
at 155 cm" 1 



TABLE I: Low-energy Raman frequencies (cm ) of 
Gd 3 N@C 2 n (40 < n < 44) below 200 cm" 1 at 90 K. Bold 
frequencies represent the Gd— cage stretching mode, which 
appears to be common for all gadolinium samples. 

Gd 3 N@C 8 o Gd 3 N@C 8 2 Gd 3 N@C 8 4 Gd 3 N@C 8 6 Gd 3 N@C 8 8 
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FIG. 3: Analysis of low-energy Raman lines of Gd 3 N@C2n 
(40 < n < 44) taken at 90 K indicating a hindered rotation 
due to the coupling of the core complex to the cage. 



When plotting this mode against the square root of the 
reciprocal reduced mass, [l/m re d U ced (GdsN — C2 n )] 1 ^ 2 , 
a linear trend is found (Fig. 4). As the carbon cage mass 
increases, a linear frequency downshift is observed due to 
a decrease in the bond strength between the core complex 
and the cage. It is interesting to note that previous X-ray 
data has distinguished GdsN^Cso, which obeys the IPR 
and contains a pyramidal GdsN unit, from GdsN@C82 
and GdsN@C84 which contain planar Gd3 units and car- 
bon cages that do not follow the IPR f7j. Although no 
X-ray or NMR data has been conducted for Gd3N@C86 
and GdaN^Css, the cages most likely follow the IPR and 
have D 3 and D 2 symmetry respectively [37], [HI . It is im- 
portant to note that unlike GdaN^Cgo, which has a pref- 
erence for icosahedral symmetry (I/,,), the other molecules 
have several IPR isomers that could be present in the Ra- 
man data. 
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FIG. 4: Analysis of Gd— cage streching mode plotted against 
the square root of the reciprocal reduced mass. As the carbon 
mass increases, there is a linear downshift in the frequency. 

Raman data of YaN^Cgo and LusN^Cgo were ana- 
lyzed in order to compare with the gadolinium complex 



Gd 3 N@C 8 o (Fig. 5). There is agreement above 220 
cm -1 , which identifies the Cso cage modes. The fre- 
quency of the M— cage stretching mode is increased to 
195.6 cm -1 for the yttrium complex, due to the change 
in reduced mass (32.7% decrease), and decreased to 155.9 
cm -1 for the lutetium complex, due to the change in re- 
duced mass (7.0% increase). Assuming harmonic motion, 
a force constant k = 12.8 ± 0.37 N/cm can be derived, 
and it can be concluded that the frequency shift of the 
M3N— Cgo oscillator is solely due to the change in metal 
mass [39j. A splitting and separation of the M— cage 
stretching modes (123.4—195.6 cm -1 ) is observed as the 
metal mass increases from yttrium to lutetium, as well 
as COM modes below 120 cm -1 . 
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FIG. 5: Comparison of Y 3 N@C 8 o, Gd 3 N@C 8 o, and 
LusN^Cso low-energy Raman data taken at 90 K. Analy- 
sis of the data identifies Cso cage modes, hindered rotation 
modes, and center of mass modes. 
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FIG. 6: Comparison of Y 3 N@C 88 and Gd 3 N@C 88 Raman 
data taken at 90 K. Analysis of the data identifies Css cage 
modes, hindered rotation modes and center of mass modes. 
Prominent peak correlations are denoted by dotted lines. 



gas, and submerged into a liquid helium storage dewar. 
The junctions were formed after thermal equilibrium was 
reached and all transport measurements were performed 
at a temperature of 4.2 K. In Fig. 7, a zero-bias conduc- 
tance peak commonly associated with Kondo scattering 
[4lj is clearly present. 

The Kondo effect is a many-body phenomenon that 
can arise from the coupling between a localized spin (An- 
derson impurity) and a sea of conduction electrons. It 
has been shown that molecules composed of magnetic 
ions can have a non-zero magnetic moment and give rise 
to Kondo scattering [II El The magnetic moment 
of the gadolinium ions is due to the localized 4/ states, 
which exhibit a non-vanishing magnetic moment and a 
localized single spin. The interaction of the / electrons 
from the three Gd 2+ ions within the Cso cage may in 
principle lead to the observed Kondo effect. 



Raman data of Y3N@Cg8 has also been studied in or- 
der to compare with GdaN^Cgs (Fig. 6). There is excel- 
lent agreement above 250 cm -1 , despite several possible 
isomers, which identifies the Css cage modes. The fre- 
quency of the M— cage stretching mode is increased to 
171.8 cm -1 for the yttrium complex, due to the change 
in reduced mass (33.3% decrease). Assuming harmonic 
motion, a force constant k = 9.9 ± 0.12 N/cm can be 
derived and one can conclude that the frequency shift 
of the M3N— Css oscillator is solely due to the change 
in metal mass [39j. Several COM modes are observed 
in Gd 3 N@C 88 at 93.5, 74.5, 60.2, 43.9, and 21.4 cm" 1 
respectively. 

B. IETS and Conductance Measurements 

The low temperature charge transport properties of 
Gd,sN@C8n .junctions were measured in a crossed- wire 
apparatus [40[. A Gd3N@Cso layer is deposited from 
a ~ 10 pM CS2 solution onto a 10 /im diameter Au wire 
and mounted into a custom-built cryostat. A second 10 
/im diameter Au wire serves as the top electrode. The 
cryostat assembly is evacuated, back filled with helium 
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FIG. 7: Experimental conductance data of Kondo effect and 
zero-bias anomaly in Gd3N@Cso taken at 4.2 K. Inset adapted 
from [4(| shows the experimental setup: Au crossed- wire ap- 
paratus forms a junction with the Gd3N@Cso thin film. 

It is important to note however, that similar Kondo 
scattering has been previously observed for junctions con- 
taining C ^p t hat do not have an intrinsic magnetic mo- 
ment |45l l46j. Therefore, it is difficult to identify if the 
spin impurity is due to the Cso cage, the gadolinium com- 
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plex, or some combination of the two. Further investiga- 
tion by magnetic STM (Scanning Tunneling Microscopy) 
or SPLEEM (Spin-polarized Low-energy Electron Mi- 
croscopy) is required to have a more detailed analysis 
of the system. 

Satellite peaks in the conductance data at ±42 mV 
(339 cm -1 ), similar to those reported [46[ for Ceo at 
±33 mV (266 cm -1 ), are also observed in the conduc- 
tance spectra of the Gd3N@C8o junction. However, un- 
like previous groups who identified the peaks as the H p (l) 
squashing mode for Ceo , the observed satellite peaks can- 
not be definitively correlated with any intramolecular 
modes of Gd 3 N@C 8 o. 
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FIG. 8: Experimental IETS spectrum obtained for 
Gd 3 N@C80 with modulation amplitude V w = 4 mV, scan rate 
= 1.5 mV/s and time constant = 1 s. The anti-symmetric 
Gd— N stretch mode is identified at 81.6 mV (658 cm -1 ) as 
well as Raman Cso cage modes at 155.1 mV (1251 cm -1 ) and 
187.5 mV (1512 cm -1 ). 

Not all junctions formed exhibited the Kondo signa- 
ture shown in Fig. 7. In fact, the transport properties 
were more or less evenly split between the Kondo reso- 
nance previously discussed and more conventional non- 
resonant tunneling and inelastic electron tunneling sig- 
natures |47j|- The inelastic electron tunneling spectrum 
of a GdaN@C8o junction is shown in Fig. 8. There is a 
strong line at 81.6 mV (658 cm -1 ) which has been identi- 
fied previously by IR spectroscopy as the anti-symmetric 
Gd— N stretching vibration of the GdsN cluster [35]. Ra- 
man Cso cage modes are also observed at 155.1 mV (1251 
cm -1 ) and 187.5 mV (1512 cm -1 ). The vibrational spec- 
trum obtained by IETS is complimentary to Raman and 
IR measurements since although there are no strong se- 
lection rules [48|, [49[ , the molecular vibrations observed 
provides insight into the conduction pathways through 
the molecular species [50|. 



IV. DISCUSSION 

In the case of Cso, its 12 pentagonal and 30 hexag- 
onal faces can be arranged in seven different structures 



consistent with the IPR [5l|. The Cgo cage has elicited 
significant interest in part because after Ceo, it has the 
smallest fuller ene cage with an icosahedral isomer. The 
isolated Cgo molecule has 234 vibrational degrees of free- 
dom: 80 x 3 = 240 degrees of freedom in total reduced 
by six corresponding to rotation and translation. Due to 
the icosahedral symmetry of the cage most of the vibra- 
tional modes are highly degenerate and fall into 62 dis- 
tinct mode frequencies. The intramolecular modes can be 
analyzed and classified according to their symmetry [52[ 
and from group theory it follows that 14 of the modes are 
Raman active (3A^ and 11H P ), 6 are IR active (Fi n ) and 
the remaining 42 are optically silent. Vibrational modes 
for the isolated 1^ — Cg^~ cage were calculated by hybrid 
DFT (B3LYP) utilizing the STO-3G/3-21G basis sets in 
the Gaussian 03 package [53[. Pure Cso cage vibrations 
should not be affected by the mass of the endohedral 
atoms [36|. These values are compared with experimen- 
tal and theoretical calculations (UFF) of Gd3N@C8o in 
Table II. 

TABLE II: Experimental and theoretical Raman frequencies 
(cm -1 ) for Gd 3 N@C 8 o, Ih-C^ , and intrinsic Gd 3 N. 

Gd 3 N@C 8 o (Exp) Gd 3 N@C 8 o (UFF) I h -C%- (DFT) 



430.4 A g (l) 428.5 A g (l) 429.7 A g (l) 
376.1 372.1 375.6 

361.3 R g (2) 370.7 H^(2) 355.3 H^(2) 

343.1 356.3 302.4 

313.2 302.2 298.1 
296.4 255.9 227.2 R g (l) 
265.7 229.1 

234.5 R g (l) 217.9 R g (l) 
165.2 Gd-cage 141.8 Gd-cage 

148.4 116.3 Gd 3 N (UFF) 

91.5 115.4 747.8 

83.6 84.2 511.5 
74.9 79.2 112.1 
59.8 75.5 94.8 
44.6 46.7 

35.8 44.3 

24.2 25.4 



The Cg^" model accounts only for the charge trans- 
fer in endohedral complexes, which is a very simplified 
approximation of the core— cage interaction. In reality, 
the electronic properties of endofullerenes are determined 
additionally by the hybridization of the core— cage or- 
bitals. Describing the presence of gadolinium only by 
adding electrons to the cage neglects some important ef- 
fects arising from Gd— cage chemical bonding. In such 
a simple model the extra electrons are delocalized over 
the whole molecule, whereas more complex calculations 
for the endofullerenes predict localization of the charge 
in the vicinity of the chemisorption site. 

The intrinsic GdaN molecule has six vibrational modes, 
two of which are degenerate. One mode with Ai sym- 
metry is a breathing mode or Gd— N symmetric stretch, 
where the Gd triangle expands and contracts. For the Ai 
wagging mode or out-of-plane breathing mode, the Gd3 
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moiety and the central nitrogen atom are displaced from 
the molecular plane in opposite directions. One of the 
E modes is a doubly-degenerate in-plane mode or asym- 
metric Gd— N stretching mode. The other E mode is 
an out-of-plane normal vibration or scissor mode, which 
involves changes in the Gd— N— Gd bond angle. All nor- 
mal mode vibrations of the free GdaN molecule are both 
Raman and IR active. 

The normal modes of the free GdaN molecule have 
been calculated by UFF. The doubly-degenerate scissor 
mode or S (Gd— N— Gd) and asymmetric stretch or v as 
(Gd-N) have values at 94.8 cm" 1 and 747.8 cm" 1 re- 
spectively. The wagging mode or 7 (Gd3— N) has the 
value of 511.5 cm" 1 , whereas the breathing mode or v s 
(Gd— N) has a value of 112.1 cm" 1 . The wagging mode 
is related to a Raman vibrational mode of GdN at 530 
cm" 1 reported previously [54j . 

V. CONCLUSIONS 

We have described the structure and vibrational spec- 
trum of the Gd 3 N@C 2n (40 < n < 44) family. The end- 
ofullerenes were characterized by Raman spectroscopy, 
IETS, UV-Vis spectrophotometry, and molecular model- 
ing (DFT and UFF). Tangential and radial cage modes, 
low-energy Gd— cage modes, intermolecular modes, and 
COM modes were identified in the Raman spectra. 
Hindered rotations were observed experimentally and 



calculated theoretically, indicating the formation of a 
GdaN— C2n (40 < n < 44) bond which reduces the sym- 
metry of the cage. As the carbon cage mass increases, 
a linear frequency downshift of the hindered rotation 
modes is observed due to a decrease in the bond strength 
between the core and the cage. Raman studies were 
conducted on Y 3 N@C 8 o, Lu 3 N@C 8 o and Y 3 N@C 8 8 for 
comparison with the Gd3N@C2 n family. GdaN^Cgo was 
further investigated by conductance measurements in a 
cross-wire apparatus at low temperature. A zero-bias 
conductance peak associated with Kondo scattering was 
observed. Based on these observations, we suggest that it 
is possible to utilize magnetic endofullerenes in magneto- 
electronic device applications, while characterizing the 
molecules through Raman and IETS. 
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